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Abstract 

We discuss the evolution of Ef small x, emphasizing the uncer¬ 
tainties related to expansion, fitting, renormalization and factorization 
scheme dependence. We find that perturbative extrapolation from the 
measured region down to smaller x and lower may become strongly 
scheme dependent. 
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It is now well established!]^ that the behaviour of F 2 in the region of small 
X and large accessed by the HERA experiments|0 provides a conhrmation 
of the double scaling behaviour |Q] predicted asymptotically in perturbative 
QCD0. As the experimental accuracy improves, it is now possible to test 
the theory beyond this simple leading order prediction, by comparing the 
data to a full next-to-leading order (NLO) determination of the x and 
dependence of F 2 . This, however, requires a study of the renormalization 
and factorization scheme dependence which characterizes perturbative com¬ 
putations, and which at small x become particularly signihcant, due to the 
growth of anomalous dimensions and coefficient functions. Moreover, the 
presence in the problem of two large scales and s = — x)/x) re¬ 

quires the choice of an expansion scheme which sums up all the appropriate 
leading (and subleading) logarithms. Here we assess the size of these am¬ 
biguities and in particular discuss how they affect the computation of F 2 
and, conversely, the extraction from F 2 of information on the form of parton 
distributions at small x. 

Determining the evolution of structure functions by solution of the renor¬ 
malization group equations in leading order corresponds to summing all logs 
of the form af(logQ^)'^(log i)'’ with p = q and 0 < r < p; double scal¬ 
ing is a consequence of the dominance at small x of the contributions with 
r = p = q, i.e., such that the two large logs are treated symmetrically. It 
is in fact possible [§] to reorganize the perturbative expansion in such a way 
that the full LO contribution to anomalous dimensions treats the two logs 
symmetrically, i.e. such that in LO each power of a* is accompanied by ei¬ 
ther of the two logs (that is, such that l<q<p, 0<r<p, 1 < p < q + r). 
This expansion scheme (the double leading scheme) can then be extended to 
NLO and beyond. Within any given scheme the structure functions are ex¬ 
pressed as power series in even though solving the renormalization group 
equations sums contributions involving large logarithms to all orders in ag. 

Consistent solution of the evolution equations in any specihed expansion 
scheme and to a given order is then (at least in principle) always possible. 
In practice, the anomalous dimensions are known through their Laurent ex¬ 
pansion in N. All the LO coefficients of this expansion are known for the 
2x2 matrix of singlet anomalous dimensions]^, but the NLO coefficients 
of the singular terms in N are only known for 7 ^ and 7 |]f 0; the corre¬ 
sponding coefficients in 7 ^ and 7 ^ can however be hxed by requiring mo¬ 
mentum conservation^]. We will henceforth consider NLO computations in 
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the double-leading scheme.^] More specifically perturbative evolution is per¬ 
formed in the usual loop expansion scheme down to a certain xq, and then 
the double-leading scheme used below it. The value of the parameter xq can 
only be determined by comparison to experiment. Here we will consider two 
extreme double-leading NLO scenarios, namely Xq smaller than any value of 
X covered by the HERA data, (i.e., in practice, Xq = 0 or ordinary two-loop 
evolution), and Xq = 0.1. 

Once Xo is fixed, renormalization and factorization schemes still have to be 
specified in order to perform NLO computations. Without loss of generality 
the renormalization scheme will be chosen to be MS: other renormalization 
schemes then correspond simply to a change of renormalization scale. The 
choice of factorization scheme is more complex. Firstly, we have a choice 
between schemes in which to all orders F 2 is directly proportional to the 
quark distribution (parton schemes, such as the DIS scheme), and schemes 
where (starting at NLO) F 2 receives a gluon contribution (such as the MS 
scheme). Furthermore, in parton schemes the coefficients of NLO singular 
contributions to the singlet anomalous dimensions also depend on the choice 
of factorization scheme: when computed in the DIS scheme 0 they contain 
a process-independent singularity in the quark sector, which is removed if 
off-shell factorization is used instead (QqDIS scheme) [|I0| . It is even possible 
to set the NLO singularities in the quark sector to zero, thereby factorizing 
the entire singularity into the starting distribution (SDIS scheme) ||1 l|| .f| Cor¬ 
responding MS schemes may be constructed by insisting that the anomalous 
dimensions be the same as those in the standard MS scheme [0, but with 
the coefficients being adjusted accordingly (so that in particular in QqMS 
scheme the process - independent singularity is removed from the coefficient 
function). 

Finally, there is still an ambiguity in the definition of the initial parton 
distributions (a ‘fitting scheme’ ambiguity), related to the fact that these can 
be fitted in a parton scheme or in an MS scheme regardless of which scheme 
is chosen to evolve. Besides providing information on the dependence of 

*Notice that this is not quite the same as the approach of ref. j^, where the higher 
order singularities are simply added to the one and two loop anomalous dimensions: in the 
double-leading expansion all the NLO terms may be treated consistently, by linearizing 
them in order to avoid spurious sub-subleading terms. 

Tt turns out that in the HERA region the results obtained in the SDIS scheme are 
essentially identical to those found by ordinary two loop evolution|^. 
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norm.(%) 

A. 

A. 



96 

103 

-0.23 

± 

0.05 

0.10 ±0.07 

57.3* 

a) 

97 

103 

-0.24 

± 

0.05 

0.12 ±0.08 

57.6 


94 

101 

-0.24 

± 

0.09 

-0.52 ±0.23 

59.2 


95 

101 

-0.25 

± 

0.10 

-0.49 ±0.26 

59.0* 


96 

102 

-0.25 

± 

0.02 

0.03 ±0.16 

64.5* 

b) 

97 

104 

-0.25 

± 

0.02 

-0.08 ±0.01 

58.1 


97 

104 

-0.12 

± 

0.02 

-0.01 ±0.20 

62.5 


97 

103 

-0.13 

± 

0.07 

-0.36 ±0.24 

57.9* 


96 

102 

-0.26 

± 

0.02 

0.12 ±0.17 

72.6* 

c) 

98 

106 

-0.24 

± 

0.02 

-0.17 ±0.09 

65.0 


97 

103 

0.10 ±0.06 

0.01 ±0.37 

73.1 


95 

101 

-0.03 

± 

0.03 

-0.75 ±0.05 

62.3* 


100 

106 

-0.13 

± 

0.04 

0.18 ±0.01 

63.4 

d) 

94 

100 

-0.22 

± 

0.06 

-0.10 ±0.10 

57.9 


100 

107 

-0.16 

± 

0.11 

0.07 ±0.25 

63.9 


89 

95 

-0.28 

± 

0.05 

-0.70 ±0.12 

73.9 


Table 1: Fitted parameters for: a) xq = 0; b) xq = 0.1, Qq factorization; 
c) Xq = 0.1, standard factorization. In each case the four entries correspond 
respectively to DIS distributions (DIS and MS evolution); MS distributions 
(DIS and MS evolution). The two entries d) show the effect on the hrst and 
fourth entry of the table of varying the renormalization scale by a factor of 
two either side. 


the results for F 2 on the specihc choice of parton parametrization, varying 
the htting scheme demonstrates the implicit scheme dependence of the htted 
parameters. 

The results of htting F 2 to HERA data Q] are summarized in the table and 
displayed in the hgure.| The free parameters are the normalizations of the 
two data sets and the small-x exponents of the quark and gluon distributions, 

Mhe corresponding results of ref. are determined by using a slightly different treat¬ 
ment of thresholds: here continuity of F2 is imposed (continuity of DIS distributions) 
whereas there continuity of the MS parton distributions was required instead. The slight 
variation of the results gives a feeling for the corresponding uncertainty. 
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which behave as as x ^ 0; the resulting (for 120 d.f.) is also given. 
All hts are performed with as{Mz) = 0.120[^; initial parton distributions 
are given at 2 GeV for the Xq = 0 hts and 3 GeV for xq = 0.1. | 

The results can be summarized as follows: 

a) Whereas the inclusion of two loop corrections improves signihcantly the 
agreement of F 2 with the data, going over to the double leading scheme has 
very little effect. 

b) Gonsequently, the data cannot yet hx the value of xq, however if xq is as 
large as 0.1 they favour Qq factorization over the standard one. 

c) In general both the relative and absolute sizes of Ag, \g depend strongly on 
expansion, htting, renormalization and factorization schemes. In particular 
if Xo = 0 in MS htting Ag ~ Ag (within errors), but in DIS htting Ag > Ag; 
while if Xo = 0.1 in MS htting Ag > Ag, but in DIS htting \ ~ Ag.0 The 
exception is that Ag when htted in DIS is independent of expansion, renor¬ 
malization and factorization scheme, since it is directly related to the small-x 
behaviour of a physical observable, F 2 . 

d) The scheme dependence of Ag is least severe in Qq factorization, where 
the gluon is generally rather hat. This provides phenomenological support 
to the theoretical expectation]^ that the input to perturbative evolution is 
of more direct physical signihcance in this scheme JJ] 

e) Whereas at hxed starting scale schemes with larger xq tend to have some¬ 
what smaller values of A (i.e. less singular inputs) the main ehect of going 
over to the double leading scheme is to reduce the sensitivity to the starting 
distribution: double scaling then results /^from a rather wide range of bound¬ 
ary conditions. 

f) Gonversely, there is a very large scheme dependence at large p (i.e. close 
to the boundary of perturbative evolution) which may signal a breakdown 
of leading-twist perturbative calculations there. This makes a perturbative 

§The starting scale should also be treated as a free parameter; it turns out however 
that a good fit can be obtained within quite a wide range of values of Qqj the resulting 
values of A being decreasing functions of Qo- 

^This seems to disagree with ref. |]^ where (on the basis of an MS calculation at two 
loops) it is claimed that Ag is significantly smaller than Xq-. it also suggests that some of 
the assumptions made in the discussion of the relative size of Ag and Ag in ref.|]lT| are 
incorrect. 

II The “initial Pomeron” reconstructed from the best-fit initial gluon distribution ac¬ 
cording to ref.]^ appears then to be soft. 
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reconstruction of the input parton distribution (and in particular the input 
gluon) from a measurement of the evolved structure function very difficult. 
Which is as it should be: evolving to smaller x and/or lower leads one 
eventually into the intrinsically nonperturbative region, 
g) Direct measurements of F 2 at larger values of p may help to reduce this 
ambiguity (or at least postpone it to yet larger values of p) by putting con¬ 
straints on xq. However, if the new data deviate strongly from the two loop 
curve this might suggest a breakdown of leading twist perturbation theory 
in this region. 

Finally we note that when the physical parameter as is also included in 
the fit, its value turns out to be largely insensitive to all of these scheme 
ambiguities, thereby allowing a determination of it from small-x structure 
function data alone [||. 

Acknowledgements: S.F. thanks G. Altarelli, S. Catani, A. Cooper-Sarkar, 
F. Hautmann, A. Martin, R. G. Roberts and A. Vogt for interesting discus¬ 
sions during the conference. 
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Figure 1: Scaling plots corresponding to double scaling (dotted); two loops 
(double leading with xq = 0) MS or DIS (solid); double leading MS (dot 
dash), QqMS (double-dot dash), QqDIS (long dashes), DIS (short dashes). 
The double scaling curve is from ref. [4], the other curves correspond to entries 
denoted by * in the table. 
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